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ABSTRACT 

We present a new determination of the mass function of galaxy clusters, based on optical virial mass 
estimates for a large sample of 152 nearby (z < 0.15) Abell-ACO clusters, as provided by Girardi et al. 
(1998). This sample includes both data from the literature and the new ENACS data. The resulting 
mass function is reliably estimated for masses larger than Mu m ~ 4 x 10 14 /i _1 M Q , while it is affected 
by sample incompleteness at smaller masses. 

We find N(> M Km ) = (6-3 ± 1.2) 10 -6 (h^Mpc)- 3 for cluster masses estimated within a 1.5 
hr x Mpc radius. Our mass function is intermediate between the two previous estimates by Bahcall 
& Cen (1993) and by Biviano et al. (1993). 

Based on the Press-Schechter approach, we use this mass function to constrain the amplitude of 
the fluctuation power spectrum at the cluster scale. After suitably convolving the PS predictions with 
observational errors on cluster masses and CO-BE-normalizing the fluctuation power spectrum, we find 
o- 8 = (0.60 ± 0.04)^0 °- 46+0 09n ° for flat low-density models and cr 8 = (0.60 ± 0.04)fV°- 48+0 - 17n ° for open 
models (at the 90% al.). 
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According to the scenario of hierarchical formation of 
G cosmic structures, clusters of galaxies arise from rare high 
Q-i peaks of the initial density fluctuation field and, at present 
time, represent the largest virialized cosmic structures. As 
a consequence, standard analytical arguments based on 
the Press & Schechter (1974; hereafter PS) approach show 
that their abundance is highly sensitive to the amplitude 
of the density fluctuations on the cluster mass scale (e.g., 
White, Efstathiou & Frenk 1993). In particular, if M is 
the typical cluster mass and Slo is the cosmic density pa- 
rameter, then the typical size of the density fluctuation 
collapsing into a cluster is R oc (fioM) 1 / 3 . Typical cluster 
masses, M ~ 5 X 10 14 /i _1 M Q , are rather close to the 
average mass contained within a sphere of 8/i _1 Mpc ra- 
dius. Therefore, it is common to express constraints from 
the observed local cluster mass function (i.e., the number 
density of clusters of a given mass), n(M), in terms of f^o 
and as, the r.m.s. density fluctuation within a sphere of 
Sh^ 1 Mpc radius (e.g., Eke, Cole, & Frenk 1996; Viana 
& Liddle 1996). 

From the theoretical side, the mass function of clusters 
is easily obtainable either from N-body simulations or, 
to a comparable degree of reliability (e.g., Lacey & Cole 
1993; Eke et al. 1996; Borgani et al. 1998), by applying 
the analytical PS recipe. On the contrary, predictions of 
cosmological models about quantities which are more di- 
rectly connected to observations, like the internal velocity 



dispersion of member galaxies, a v , and the AT -ray temper- 
ature of the intracluster gas, Tx, are much less straight- 
forward and rely on assumptions which are quite debated. 
As for a v , N-body simulations allow reliable predictions 
of the internal velocity dispersion of clusters, but only for 
the DM particles (e.g., Frenk et al. 1990; Borgani ct al. 
1997; van Haarlem, Frenk, & White 1997), while the con- 
nection to the observed a v requires a plausible model of 
galaxy formation. Numerical simulations by Frenk et al. 
(1996) have shown that present uncertainties about the 
physics of galaxy formation in clusters can lead to an in- 
correct determination of the virial mass by up to a factor 
five. As for Tx, it is generally believed to be a rather 
accurate indicator of the cluster mass (e.g., Evrard, Met- 
zler, & Navarro 1996). On the other hand, both analytical 
predictions and hydrodynamical cluster simulations, based 
only on the adiabatic physics of the intracluster gas, are 
able to reproduce neither the observed relation between 
bolometric A" -ray luminosity and temperature (e.g., Eke, 
Navarro & Frenk 1997; Bryan & Norman 1998, and refer- 
ences therein), nor the observed size of core radii in the 
gas distribution (e.g., Makino, Sasaki, & Suto 1998, and 
references therein). Such a discrepancy creates a need to 
introduce additional physics to explain the thermodynam- 
ics of the the intra-cluster gas (e.g., Bower 1997; Cavaliere, 
Menci, & Tozzi 1997, and references therein). 

From the observational side, cluster masses are inferred 
from either X-ray or optical data, under the general hy- 
pothesis of dynamical equilibrium. Estimates based on 



1 Here h is the Hubble constant in units of 100 km s x Mpc 
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gravitational lcnsing do not require assumptions about the 
dynamical status of the cluster, but a good knowledge 
of the geometry of the potential well is necessary (e.g., 
Narayan & Bartelmann 1997). Claims for a discrepancy 
(by a factor 2-3) between cluster masses obtained with dif- 
ferent methods casted doubts about the general reliability 
of mass estimates (e.g., Wu & Fang 1997). However, re- 
cent analyses have shown that such discrepancies can be 
explained by the different way in which strong cluster sub- 
structures bias mass estimates based on different methods 
(Allen 1997; Girardi et al. 1997). In particular, the anal- 
ysis by Girardi et al. (1998, hereafter G98) demonstrated 
that, for nearby clusters without strong substructures, a 
good overall agreement exists between X-ray and optical 
virial mass estimates. 

A reliable determination of cluster masses is a necessary 
but not sufficient condition to guarantee a reliable esti- 
mate of the mass function. Indeed, the latter also requires 
a cluster sample which is large enough to be a fair repre- 
sentation of the cluster population in the local Universe, 
at least above a certain mass limit. In past years, such 
requirements limited the possibility of obtaining a stable 
result about the cluster mass function and, in fact, the 
only two estimates so far presented (Bahcall & Ccn 1993, 
hereafter BC93; Biviano et al. 1993, hereafter B93) turn 
out to be quite discrepant. BC93 derived cluster masses 
by using global scaling relations to connect mass to clus- 
ter richness and X-ray temperature. B93 estimated clus- 
ter masses from a homogeneous virial analysis of about 70 
clusters and found a larger abundance, by about one order 
of magnitude, at the mass scale M~5x 10 14 /i _1 M Q (cf. 
also Figure 3 below). 

A new and more robust determination of the observa- 
tional mass function based on optical virial analysis is 
now possible thanks to the recent availability of a large 
amount of rcdshift data for rich clusters, as provided by 
the ESO Nearby Abell Cluster Survey (ENACS; Katgert 
et al. 1996; 1998). Recently, G98 computed virial masses 
in a homogeneous way for a large sample of ~ 170 clusters, 
obtained by combining the ENACS sample with data from 
other authors (see also Fadda et al. 1996; hereafter F96). 
The extension of this sample makes it the largest available 
data set over which to reliably estimate the cluster mass 
function. 

In the following analysis we will use the G98 results with 
the twofold aim of (a) estimating the optical virial mass 
function of nearby galaxy clusters, and (b) investigating 
the subsequent constraints on the power spectrum of cos- 
mic density fluctuations. 

This paper is organized as follows. In § 2 we briefly de- 
scribe the statistical sample of clusters masses as obtained 
by G98 and then present the determination of the cluster 
mass function. Afterwards, we compare it with previous 
determinations based on optical data by BC93 and B93. 
In § 3 we apply the PS approach to place constraints on 
the fluctuation power spectrum for different cosmological 
models. We give a brief summary of our main results and 
draw our conclusions in § 4. 



2 THE DETERMINATION OF THE MASS 
FUNCTION 

2.1 The Statistical Sample of Cluster Masses 

In the following analysis we largely use the results by 
G98, who analyzed optical data (galaxy positions and ve- 
locities) for 170 nearby clusters having z < 0.15. The 
original optical data partly comes from the literature and 
partly from the ENACS dataset (Katgert et al. 1996; 
1998). In this section we summarize the G98 results which 
are relevant for the mass function analysis and we refer to 
that paper for further details. 

The total mass of each cluster is computed by apply- 
ing the virial theorem to the member galaxies, under 
the usual assumption that the mass distribution follows 
the galaxy distribution (e.g., Limber & Mathews 1960; 
Giuricin, Mardirossian, & Mezzetti 1982; Merritt 1988; 
B93). This assumption is supported by several pieces of 
observational evidence, coming both from optical (e.g., 
Carlberg, Yee, & Ellingson 1997) and X-ray data (e.g., 
Watt et al. 1992; Durret et al. 1994; Cirimele, Nesci, & 
Trcvcsc 1997), as well as from gravitational lensing data 
(e.g., Narayan & Bartelmann 1997). Although some level 
of uncertainty still remains about the very central cluster 
regions, detailed optical analyses (Cirimele et al. 1997; see 
also Merritt & Tremblay 1994 for the Coma cluster) sug- 
gest a very peaked galaxy distribution, in agreement with 
data on the dark matter distribution from gravitational 
lensing (see, e.g., Narayan & Bartelmann 1997). 

In order to fit the galaxy distribution for each cluster, 
G98 used a King-like profile, p(r) oc (1 + r 2 /rl)~ 3 ^/ 2 , 
where r c is the core radius and [3 is a free exponent. They 
found that, on average, p(r) oc r~ 2A for the scaling of the 
galaxy number density in the outer part of the clusters, 
thus in agreement with previous studies (e.g., Bahcall & 
Lubin 1994; Girardi et al. 1995), with a peaked distri- 
bution in the inner part (r c ~ 50/i _1 kpc; see also § 8 of 
G98). 

G98 computed cluster masses within the radius 0.002 -a v 
h^ 1 Mpc, that represents a rough estimate of the virial- 
ization radius R vir , within which dynamical equilibrium is 
expected to hold. A surface term correction has also been 
applied to the usual virial mass estimator which, other- 
wise, would have overestimated the true cluster mass (e.g., 
The & White 1976; Binncy & Trcmainc 1987; Carlberg et 
al. 1997b). This correction amounts to 19%, very similar 
to that suggested by Carlberg et al. (1997b) for CNOC 
clusters. 

G98 compared their optical mass estimates to those de- 
rived from A -ray analyses for a list of 66 clusters, which 
they compiled from the literature, and found an overall 
good agreement (cf. Figure 5 of G98). This agreement 
is expected in the framework of the two assumptions that 
mass follows the galaxy distribution and that clusters are 
not far from dynamical equilibrium (see also Evrard et al. 
1996; Schindler 1996; Carlberg et al. 1997c). Only about 
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10% of clusters show at least two strongly superimposed 
peaks in their velocity distribution thus indicating that 
they probably are far from dynamical equilibrium. The 
ambiguity in the treatment of these clusters - peaks can 
be treated as either disjoined or joined into a single struc- 
ture - leads to mass estimates which differ by a factor ~ 3 
on average (cf. § 6 of G98). The effect of these clusters 
with uncertain dynamics in the estimate of the mass func- 
tion will be discussed in the following. 

As a concluding remark, we point out that the size of our 
sample, along with the overall good agreement between X- 
ray and optical cluster masses, makes it the largest avail- 
able data set over which to reliably estimate the cluster 
mass function. 

2.2 From Masses to the Mass Function 

The results from G98 that we will use in our analysis of 
the cluster mass function are the following: (i) the values 
of the (line-of-sight) velocity dispersions u v {op in G98); 
(ii) the values of the "corrected" virial masses, M {Mqv hi 
G98), contained within the radius 0.002 -a v h~ x Mpc; (Hi) 
the King-like galaxy number-density profile which is fitted 
for each cluster and used to rescale M at different radii. 
Hereafter, we limit our analysis to the 152 Abell-ACO 
clusters which span the wide interval of richness classes 
— 1 < R < 4 (cf. Tab. 1 of G98), where the richness class 
R — — 1 is assigned to clusters with Abell number counts 
N c < 30 (Abell, Corwin, & Olowin 1989). 

The G98 sample is complete neither in volume nor in 
richness. However, there exists a well defined, although 
rather scattered, correlation between cluster richness and 
mass. Indeed, in Table 1 we show that, although the clus- 
ter mass increases with richness, there is a large mass su- 
perposition between different richness classes R at a given 
mass. This is also explicitly shown in Figure 1, where we 
plot the mass distribution for different R classes. Mazure 
et al. (1996, § 6; M96 hereafter) showed that the remark- 
able scatter in the a v -R relation is largely intrinsic. The 
propagation of this scatter is the main reason for the scat- 
ter in the M-R relation since the relation between mass 
and a v is rather strict (with a scatter of ~10%; see the 
upper panel of Figure 2). Therefore, we expect that the 
scatter in the M-R relation is also largely intrinsic, much 
like that in the a v -R relation. 

TABLE 1 

M vs. Abell Richness R 



R 


Cluster # 


log(M/M @ ) 


-1 


7 


14.08 ±0.15 





20 


14.13 ±0.10 


1 


74 


14.24 ±0.06 


2 


35 


14.61 ±0.08 


>3 


7 


14.74 ±0.11 



The 17 multipcaked clusters (see text) are not reported here. 

In order to obtain a more representative cluster sam- 
ple which also accounts for the presence of this scatter in 



the M-R relation, we resample our clusters so as to mimic 
the i?-distribution of the Edinburgh-Durham Cluster Cat- 
alogue (EDCC), which is claimed to be complete also for 
rather poor clusters (Lumsden et al. 1992). We compute 
the mass distribution by using 10,000 random extraction 
of the measured M values, distributed according to the 
EDCC ^-distribution, by rescaling the EDCC richness to 
Abell's (see also B93; Girardi et al. 1993; F96). The effec- 
tiveness of this resampling procedure has been shown by 
F96, who applied it to a sample similar to the G98 one in 
order to obtain the cr„-distribution. As a result (cf. their 
Figure 6), they were able to reproduce the a v -distribution 
derived by M96 for the volume- and richness-complete 
ENACS sample, with the further advantage of extending 
the richness range from R = 1 to R = — 1 and, therefore, 
the a v range of completeness from u v > 800kms _1 down 
to u v > 650 km s -1 . 
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Fig. 1. — The mass distributions for clusters belonging to dif- 
ferent richness classes. R = and R = — 1 classes are here 
considered together, owing to their similar mass distribution 
and to the small number of R = — 1 clusters (see also text and 
Table 1). 

In the present analysis, we apply the resampling proce- 
dure both to the whole sample of 152 R > — 1 clusters and 
to the sample of 120 R > 1 clusters and obtain the corre- 
sponding mass distributions. Since no substantial differ- 
ence exists between masses for R = — 1 and R = classes 
and owing to the small number of R = —1 clusters, we 
treat these two classes together. We verified that no dif- 
ference is found in the final results if the two classes are 
instead treated separately. 

Having determined the shape of the mass function 
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through this resampling procedure, we have still to fix 
its amplitude by resorting to an external normalization. 
To this purpose, we adopt the cluster volume density 
N = 8.6 x 10- 6 (h- 1 Mpc)- 3 for R > 1 clusters, as 
provided by M96 for ENACS, scaled to the EDCC R- 
frequencies. Since the M96 value was corrected for the 
incompleteness of the Abell-ACO catalog with respect to 
the EDCC catalog, our normalization is consistent with 
the resampling procedure. 




10- 8 
10 13 



10 14 
M (h" 1 M ) 



Fig. 2. — The lower panel shows the cumulative mass func- 
tion for R > —1 clusters (filled points with solid line) and for 
R > 1 clusters (open points with dashed line). Error bars 
represent 1-a uncertainties. The two curves above and below 
the R > — 1 mass function correspond to two different ways 
in which clusters, which appear as multi-peaked in velocity 
space, are treated (see text) . The mass completeness limits for 
the two richness classes (indicated by the vertical lines) are in- 
ferred from the relation between mass and velocity dispersion, 
as reported in the upper panel and from the corresponding a v 
completeness limits. 

We estimate the mass-completeness limit in our sam- 
ple from that in u v . The broad relation between richness 
and a v limits the completeness in velocity dispersion to 
a v > 650 km s -1 for a sample with R > — 1 clusters (F96) 
and to a v > 800 km s -1 for the sample with R > 1 clus- 
ters (M96). Then, the less scattered M-a v relation (see 
the upper panel of Figure 2) is used to estimate the mass 
completeness limit. Wc find M> 4 x 10 14 h^ 1 M Q for 
R > -1 and M> 7 x 10 14 /i" 1 M & for R>1. 

In the lower panel of Figure 2 we show the resulting 
cumulative mass functions, N(> M) = n(M)dM for 



R > — 1 clusters (filled points) and that for R > 1 clusters 
(open points). As expected, the curves flatten at small 
masses due to the effect of the incompleteness, while they 
are almost overlapping within the common completeness 
range. Errors in the mass function have been estimated by 
taking the contribution from (a) the Poissonian uncertain- 
ties from the finite statistics in the G98 sample; (b ) the un- 
certainty in the resampling procedure connected with the 
Poissonian errors in the richness frequency distribution of 
the EDCC; (c) the uncertainties in the individual cluster 
mass estimates. In this plot and in the following, errorbars 
for the mass function correspond to the 1-a overall errors 
contributed by such three sources. 

The ambiguity in how to treat the 17 highly substruc- 
tured clusters found in the sample give rise to two dif- 
ferent determinations of the mass function, depending on 
whether each of them is treated as a single structure or as 
disjoined into two objects (cf. also § 6 in G98). The re- 
sulting mass functions for the R > — 1 case are plotted in 
Figure 2 with the lower and upper light continuous lines, 
respectively. It turns out that they are always very close, 
the difference being within the statistical uncertainties. In 
the following we will refer to the N(> M) plotted with 
the heavy solid line, which is computed by assigning equal 
weights to joined and disjoined peaks. 

2.3 Comparison with Previous Estimates 

In order to compare our optical virial mass function with 
those provided by BC93 and B93, we rescale our cluster 
mass estimates to the same 1.5 h" 1 Mpc radius by using 
the fitted galaxy distributions (see the discussion in § 2.1) 
and recompute the mass function. The result of this com- 
parison, which is plotted in Figure 3, shows that our mass 
function is intermediate between the two previous deter- 
minations. Following a maximum-likelihood approach, we 
find the following expression for our mass function: 

n(M)=n* f^-V'e-M/M* (1) 



with n* = 2.6; 
M* = 2.6+ ' 8 



x 10 



-0.5 
-0.4 

x 10 14 /i- 



^(h^Mpc)- 



^IO^H^Mq)- 1 and 



J -0.6 ~ iu " lrl Q- 
The mass function by B93 (upper panel of Figure 3) is 

based on cluster masses estimated from a homogeneous 
virial analysis and on the same resampling procedure that 
we followed here. Therefore, it is directly comparable to 
our results. Since B93 did not provide a normalization for 
their mass distribution, we fix it according to the R > 1 
cluster number density provided by M96. The discrepancy 
with respect to B93 can be explained by the fact that the 
G98 mass estimates are on average 40% smaller than their 
estimates. This is explicitly shown by the dashed curve 
in the plot, which gives the mass function of B93 after 
just translating masses by 40% toward smaller values. As 
discussed in detail by G98 (cf. their § 8), two factors 
contribute to this difference. First, the more rigorous al- 
gorithm for removing interlopers gives a smaller o~ v by a 
factor of 10%, thus translating into a 20% difference in 
mass. Second, the use of the surface correction term in 
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the virial theorem (cf. § 2.1 and Carlberg et al. 1997b) 
decreases the mass estimate by about a further 20%. 
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Fig. 3. — The comparison between this and previous determi- 
nations of N(> M) by Biviano et al. (1993; upper panel) and 
by Bahcall & Cen (1993; lower panel; only their determination 
from optical cluster data is reported here). The vertical line in 
both panels indicates our limiting completeness mass. In the 
upper panel the dashed line indicates the B93 mass function 
after rescaling their masses by 40% to smaller values. In the 
lower panel the four triangles joined by a dashed line indicate 
the determination of iV(> M) we would have obtained from 
our sample by neglecting the scatter in the mass-richness re- 
lation (see text). The plotted errorbars for iV(> M) represent 
l-cr uncertainties. 

As for the optical data analysis by BC93, they associ- 
ated the Abell or EDCC cluster richness to mass according 
to a suitable one-to-one scaling relation (cf. their eq. 4; 
note that this relation falls within the range, M(< 1.5 
h- 1 Mpc) = (0.5 - 0.7) 10 13 A C , that we find to hold for 
our cluster sample). The resulting N(> M) turns out 
to be somewhat smaller than ours, the difference being 
larger at larger masses (cf. the lower panel of Figure 3). 
The main reason for this difference is due to the fact that 
the procedure by BC93 amounts to neglect any intrinsic 
scatter in the mass-richness relation (cf. § 2.2 about the 
necessity of accounting for this scatter). In order to show 
the effect of neglecting this scatter on the BC93 results, 
we consider the mass limits corresponding to the R > 0, 
R > 1, R > 2 and R > 3 clusters, according to the BC93 
mass-richness relation. At these mass limits, we show the 
fractions of our resampled clusters (triangles connected by 



a dashed line in the lower panel of Figure 3) corresponding 
to the four above richness ranges, when neglecting the con- 
tribution from poorer clusters. This shows the relevance 
of properly including the effect of the scatter in the M-R 
relation and by how much the mass function would have 
been underestimated without accounting for it. 

In this context, the normalization is only a minor source 
of discrepancy between BC93 and our N(> M) since 
they are comparable at very small masses (~ 1 — 2 x 
10 /i _1 Mq) when the whole cluster sample is considered. 
However, we stress that, at the small mass values which 
are beyond our estimate of completeness limit, both BC93 
and our N(> M) neglect the significant contribution from 
the poor galaxy systems which are not recovered in Abell 
and EDCC catalogues. 

3 IMPLICATIONS FOR COSMOLOGICAL 
MODELS 

In order to compare our mass function to predictions 
of models for cosmic structure formation, we follow the 
Press-Schechter (1974, PS) approach, which provides the 
following expression for the number density of local clus- 
ters with mass in the range [M, M + dM] : 



n{M) dM 



2 p_ 

7T M 2 



CM 



dlog a M 



d\ogM 



exp 



2^m 



dM. 

(2) 

Here p is the average matter density and S c is the linear- 
theory overdensity for a uniform spherical fluctuation 
which is now collapsing; it is S c = 1.686 for f2o = 1, 
with a weak dependence on f2o for both flat and open 
geometries (e.g., Eke et al. 1996). The r.m.s. lin- 
ear density fluctuation om at the mass scale M is re- 
lated to the fluctuation power spectrum according to 
a 2 M = (27T 2 )- 1 J^dkk 2 P(k)W 2 (kR). Here W(x) = 
3 (sin x — x cos x) /x 3 is the Fourier transform of the window 
function, assumed to have the top-hat profile. Accord- 
ingly, the mass M of a cluster arising from the collapse of 
a fluctuation with typical size R is M = (4tt /3)pR 3 . The 
PS mass function has been compared with N-body simu- 
lations by several authors (e.g., Eke et al. 1996; Lacey & 
Cole 1996; Borgani et al. 1998; Gross et al. 1998) and has 
been generally shown to provide a rather accurate descrip- 
tion of the abundance of virialized halos of cluster size. 

Herebelow, we assume for the power spectrum of den- 
sity fluctuations the expression P(k) = AkT 2 (k), where 
the transfer function is given by 



T{q) = 



ln(l + 2.34g) 



[1 + 3.89g 



2.34g 
(16.1g) 2 + (5.46g) 3 ^ 



-1/4 



(3) 



with q — k/Th. This expression has been provided by 
Bardeen et al. (1986) for CDM models with negligible 
baryon contribution, once T = Q h is taken for the shape 
parameter. In the following, however, we will take T to 
be a free quantity to be fitted to observations, thus inter- 
preting eq. (^) as a purely phenomenological expression for 
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the transfer function. As for the amplitude of the power- 
spectrum, it is customary to express it in terms of as, the 
r.m.s. fluctuation amplitude within a top-hat sphere of 
8/i _1 Mpc radius. Therefore, the class of models that we 
are considering is represented by three parameters, CIq, T 
and as, either with or without a cosmological constant 
term, J7a = 1 — to provide flat spatial geometry. If 
we further require our model to satisfy the constraint pro- 
vided by the four-year COBE data (see, e.g., Bunn & 
White 1997 and Hu & White 1997, for flat and open ge- 
ometry, respectively), then we would have a one-to-one 
relation between the shape T and the amplitude as for 
each value of Oo, thus reducing to two the number of free 
parameters. 

The cluster mass which appears in the PS formula refers 
to the mass M„; r contained within the virialization radius 
R v ir, whose value depends on the cosmological parameters 
Qq and CI a- Therefore, for the sake of comparison with PS 
predictions, the mass function presented in the § 2.2 should 
be suitably re-estimated model by model. If A c is the ra- 
tio between the mean density within the virial radius and 
the critical average density, whose value depends on Slo 
and fl A (e.g., Kitayama & Suto 1996; A c = 18tt 2 ~ 178 
for n = 1), then R 3 mr = 3tt M vir /(4TrpA c ). Therefore, 
the value of R V i r depends on M V i r , so that we resort to 
the following recursive procedure to compute it for each 
cluster. After taking the mass contained within the ra- 
dius 0.002cr„ as a first guess for M mr (cf. § 2.1), we (i) 
provide a new estimate of R V i r according to the above ex- 
pression, and (ii) estimate the new M V i r by rescaling the 
previous value to the updated R V i r , using the appropriate 
cluster density profile (cf. § 2.1). We explicitly compute 
the observational mass function n(M) for Slo = 0.2, 0.4, 1 
for both open and flat models, while it is estimated by 
linear interpolation for intermediate Qo values. In our 
analysis, we applied eq.(^) at z = 0.05, which represents 
the median redshift of the G98 cluster sample. In any 
case we find that, for the considered range of models, the 
mass function can always be well fitted by the parame- 
teric expression of eq. ([!]), with M* ~ 3.5 x 10 14 /i _1 M o 
and normalization varying in a rather small range, n* ~ 
(0.9- 1.3) x lO" 5 ^" 1 Mpc)- 3 (I0 li h- 1 M Q )- 1 . 

A further reason of caution when comparing data and 
model predictions on n(M) is related to observational un- 
certainties in estimates of cluster masses, which are in gen- 
eral not negligible. G98 provides mass uncertainties by 
propagating statistical errors in the estimates of the ve- 
locity dispersion a v and of the virial radius (cf. col. 3. of 
Table 3 in G98). Such errors are correlated with cluster 
masses, being relatively smaller for richer, better sampled 
clusters. They range from ~ 15% at M ~ 2 x lO 1B fc _1 M 
to ~ 35% at M ~ 10 u h~ 1 M e . Therefore, the measured 
mass function results from the convolution of the intrinsic 
n(M) with mass errors. In order to account for them, we 
convolve the PS mass function of eq.(^) with the uncer- 
tainties in the mass estimates, as provided by G98. 

We show in Figure 4 the observational n{M), for Qq = 
0.4 and CI a — 0.6, in the mass range where the sample is 



considered as complete (cf. § 2.2). Also plotted are the 
predictions for models with CIq = 0.4 and Cl\ = 0.6, in or- 
der to show the effect of varying as and T separately (lower 
and upper panel, respectively). From the upper panel we 
show that n(M) is rather insensitive to the shape of P(k) 
(see, e.g., White et al. 1993), at least for M< 10 15 h- 1 M e , 
where n(M) is better determined. By varying T over a 
rather large range, much larger than that allowed by the 
power-spectrum of the observed galaxy distribution (e.g., 
Peacock & Dodds 1994), does not significantly change the 
goodness of the fit. On the contrary, our mass function 
represents a strong constraint for the amplitude of the 
power-spectrum (cf. the lower panel). The heavy and 
light curves in the lower panel show the PS predictions 
with and without convolving with mass uncertainties, re- 
spectively. Although the effect of such errors is never large, 
it turns out that at the smallest mass scale its size is of 
the same order of the error bar in the observational n(M). 
Therefore, we prefer not to neglect it in our analysis. 




10 15 
M (h"'M s ) 

Fig. 4. — The effect of varying the power-spectrum shape pa- 
rameter r (upper panel) and as (lower panel) on the error- 
convolved Press-Schechter mass function. Heavy and light 
curves in the lower panel are the PS mass function with and 
without convolving the model predictions with the uncertain- 
ties in the mass estimates. Filled dots represent the obser- 
vational n(M), with masses estimated within the virialization 
radius for a model with fio = 0.4 and J1a = 0.6. 

If we impose the COBB-normalization, our observa- 
tional n(M) provides a relation between as and CIq in 
the model parameter space. In order to estimate the 
confidence level for model rejection, we adopt a x 2 - 
minimization procedure with respect to as (or, equiva- 
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lently, with respect to T) for each value of f^o- The proba- 
bility for a model to be accepted is then computed as the 
probability that the observed n(M) comes from the par- 
ent model distribution with Gaussian random variations 
in logarithmic units given by the size of the error bars. 
As a result, we find that the 90% c.l. constraints on the 
os-f^o plane can be represented by the fitting expressions 

a s = (0.60 ± 0.04) x Q-°- 4 6+o.09fi . j^a = 1 — f^o 

(7 g = (0.60 ± 0.04) X f7-0-48+0.17f2o . Q A = Q (4) 

which are accurate within 2% in the range 0.2 < Oo < 1. 
We note that these constraints come from the mass func- 
tion estimated above the mass completeness limit which 
is appropriate for R > — 1 clusters (cf. § 2.2). If, in- 
stead, we use the R > 1 mass completeness limit, Mu m ~ 
7x lO^/i^M©, then we find a s = (0.62±0.08) for fl = 1, 
with similar rescalings as in eq.(^) for Qq < 1. Our re- 
sults are therefore unaffected by narrowing the complete- 
ness mass range. The only effect is that of increasing the 
uncertainties, due to the more limited mass range over 
which theoretical predictions are tested. 

Our results can be compared with similar results ob- 
tained by other authors from the cluster abundance. 
White et al. (1993) obtained cr 8 = 0.57 ± 0.05 for fl = 1 
from the median velocity dispersion of Abell clusters, as 
provided by Girardi et al. (1993), and from the temper- 
ature functions by Henry & Arnaud (1991) and Edge et 
al. (1990). Viana & Liddle (1996) obtained a 8 ~ 0.6 
for f^o = 1 by fitting the cluster temperature function by 
Henry & Arnaud (1991) at 7 keV. A somewhat smaller 
normalization, as — 0.52 ± 0.04 for Slo = 1, has been 
found by Eke et al. (1996), also based on the analysis 
of the cluster temperatures by Henry & Arnaud (1991) 
(see also Markevitch 1998, for a similar result based on re- 
cent ASCA temperature data). Oukbir, Bartlett & Blan- 
chard (1997) found the same result as White et al. (1993), 
also from the A" -ray temperature function. Borgani et al. 
(1998) obtained a 8 = 0.58 ± 0.02 for fl = 1 from the 
local A -ray luminosity functions provided by Ebeling et 
al. (1997) and Rosati et al. (1998). It is remarkable that 
all such results, although based on largely different data 
sets are rather consistent with each other and converges 
to indicate that as — 0.6 for £Iq = 1, with rescalings for 
different Qq values which are in general quite close to those 
in eq.(|). 

4 CONCLUSIONS 

In this paper we have presented a new determination of 
the mass function n(M) of nearby galaxy clusters. This 
analysis is based on the virial mass estimates by Girardi 
et al. (1998, G98), for a sample of 152 Abell-ACO clusters 
which includes both data from the literature and the new 
ENACS data (Katgert et al. 1998). The extension of this 
sample, along with the good agreement found by G98 be- 



tween optical and virial masses, makes it the largest data 
set over which to reliably estimate the local (z ~ 0.05) 
cluster mass function. 

The main results of our analysis can be summarized as 
follows. 

(a) After applying a suitable resampling procedure to 

account for the lack of volume- and richness- 
completeness of the G98 sample, n(M) is reliably 
computed for masses larger than Mu m ~ 4 x 
10 14 /i _1 M Q . At this mass limit, the value of the 
cumulative mass function is A^(> M) = (6.3 ± 1.2) x 
10 -6 (/i -1 -/V/pc)~ 3 for masses estimated within the 
1.5ft. -1 Mpc radius. The corresponding differential 
mass function is fitted with a Schechter-like func- 
tion [cf. eq. ([l|)], with normalization n* — 2.61q4 x 
IQ-^h^Mpc)- 3 ^ 14 ^ 1 ^^)- 1 and characteristic 
mass M* = 2.6±^ x IO^/i^Mq. 

(b) A comparison with previous mass function estimates, 

also based on optical data, by Biviano et al. (1993, 
B93) and Bahcall & Cen (1993, BC93), shows that 
our n(M) is intermediate between these two, B93 
and BC93 providing a larger and smaller estimate, 
respectively (cf. Figure 3). 

(c) A comparison with the predictions from models of cos- 

mic structure formation demonstrates that our ob- 
servational mass function provides a robust deter- 
mination of the amplitude of the fluctuation power 
spectrum at the cluster mass scales (cf. Figure 4). 
After suitably convolving the prediction from the 
Press-Schechter (1974) mass function with the ob- 
servational errors in the mass estimates, we deter- 
mine as for COBEI-normalized models with a ~ 7% 
uncertainty (at the 90% c.l.). It turns out that 
a 8 — 0.6 for Q = 1, with appropriate rescalings 
at different Qq values [cf. eqs.(|])]. 

As a concluding remark, we point out that a precise 
determination of n(M) for local clusters, combined with 
a similar virial analysis of clusters at higher redshift, 
z ~ 0.3-0.5, would provide information about the degree 
of evolution of the cluster abundance and, therefore, would 
allow one to estimate as and Slo separately (see, e.g., Bah- 
call, Fan, & Cen 1997; Carlberg et al. 1997a). The im- 
plication of the mass function presented in this paper on 
the evolution of the cluster abundance will be discussed in 
detail in a future paper. 
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